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Abstract. In attempt to explore half-metallic properties of the double perovskites 
Sr2_xY x VMoOe and S^-^Y^VTcOg, we construct an effective low-energy model, 
which describes the behavior of the t2 g -states of these compounds. All parameters 
of such model are derived rigorously on the basis of first-principles electronic structure 
calculations. In order to solve this model we employ the optimized effective potential 
method and treat the correlation interactions in the random phase approximation. 
Although correlation interactions considerably reduce the intraatomic exchange 
splitting in comparison with the Hartree-Fock method, this splitting still substantially 
exceeds the typical values obtained in the local-spin-density approximation (LSDA), 
which alters many predictions based on the LSDA. Our main results are summarized 
as follows: (i) all ferromagnetic states are expected to be half-metallic. However, their 
energies are generally higher than those of the ferrimagnetic ordering between V- and 
Mo/Tc-sites (except S^VMoOg); (ii) all ferrimagnetic states are metallic (except fully 
insulating Y2VTCO6) and no half-metallic antiferromagnetism has been found; (iii) 
moreover, many of the ferrimagnetic structures appear to be unstable with respect to 
the spin-spiral alignment. Thus, the true magnetic ground state of the most of these 
systems is expected to be more complex. In addition, we discuss several methodological 
issues related to the nonuniqueness of the effective potential for the magnetic half- 
metallic and insulating states. 
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1. Introduction 

Recently double perovskites A2- x A' x BB'Oq, where A and A' are the di- and three- 
valent (typically alkali-earth or rare-earth) elements, respectively, and B (B') are 
the transition- metal elements, have attracted considerable attention pQ, following 
experimental discoveries of the room-temperature ferromagnetism and the colossal- 
magnetoresistance effect in Sr 2 FeMoOg [2], which make these materials promising in 
the field of spin electronics. The extraordinary properties of the double perovskites are 
frequently attributed to their half-metallic electronic structure, although the situation 
can be rather subtle [3]. Nevertheless, the first discoveries spark the new research 
activity aiming at the systematic search of half-metallic double-perovskite materials. 
Many such predictions has been made on the basis of first-principles electronic structure 
calculations. Besides ferromagnetic (FM) compounds, a particular attention was paid 
to the half-metallic antiferromagnets - the systems, where the half-metallicity coexists 
with the zero net magnetization [H [5]. However, even apart from the technological 
issues, related to the defects in the double-perovskite structure, there is a number of 
fundamental problems, which call many predictions into question: 

(i) Most of the calculations are based on the local-spin-density approximation (LSDA) 
or its refinement - the generalized gradient approximation, which have many 
limitations for the transition-metal oxides. Particularly, all these theories employ 
the functional dependence of the exchange-correlation energy obtained in the limit 
of homogeneous electron gas. However, many unique properties of the double 
perovskites are related to the alternation of two different transition-metal sites 
B and B', and in this sense the systems are essentially inhomogeneous. Thus, even 
though the double-perovskite compounds may be classified as itinerant electron 
magnets, it does not necessary mean that their properties are well described within 
LSDA; 

(ii) In the vast majority of calculations, it is believed that the LSDA+?7 technique, 
which incorporates the physics of on-site Coulomb interactions [HI [7], can serve as 
a reasonable alternative to LSDA and provide more superior description for the 
properties of double perovskites. However, the LSDA+t/ scheme for the double 
perovskites is ill-defined: there are many uncertainties related to the values of the 
on-site Coulomb repulsion U at the sites B and B', which are frequently taken 
as adjustable parameters, the subspace of "localized orbitals", and the double- 
counting term [§]. Moreover, the total energy in LSDA+£7 is taken in the ad hoc 
form, which is valid only in the atomic limit and equivalent to some constrained 
density-functional theory [7J . In such a situation, all predictions of the half-metallic 
behavior for the double perovskites should be taken very cautiously. 

In this paper we will illustrate how these problem could be resolved by applying 
the optimized effective potential (OEP) method, and what the consequences on the 
electronic and magnetic properties of the double perovskites could be. The OEP method 
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is formulated in the spirit of Kohn-Sham density functional theory [9], by introducing 
an auxiliary one-electron problem, which is specified by some local effective potential 
and which is used for the search of the true spin density through the minimization of the 
total energy of the system [TUJ [HI [T2J [13] . In this work we will apply the OEP technique 
for the solution of the low-energy model, where all the parameters are derived rigorously, 
on the basis of first-principles electronic structure calculations. Besides simplicity, such 
a model analysis typically allows one to concentrate directly on the microscopic picture, 
underlying the considered phenomena. For the first applications of the OEP model, it is 
more convenient to consider the t 2g -systems, whose electronic and magnetic properties 
are predetermined by the behavior of the ^-states of the sites B and B' located near 
the Fermi level. In this case, the effective potential is specified by only three parameters, 
which should be used for the minimization of the total energy. 

We select two prototypical examples Sr 2 _ x Y :r VMo06 and Sr 2 -xY x VTc06, for which 
we adopt the undistorted double perovskite structure with the lattice parameters a = 
7.846 and 7.899 A, respectively. In fact, a = 7.846 A is the experimental lattice 
constant for SrLaVMoOg [H] , which is close to theoretical values obtained in the first- 
principles calculations [15], and a = 7.899 A is the theoretical lattice constant for 
La 2 VTc0 6 [16J. Deviations from the ideal double perovskite structure are expected to 
be small [16] and were not considered here. The band-filling dependence due to the 
substitution of Sr by Y was treated in the virtual-crystal approximation. According to 
the electronic structure calculations in the local-density approximation (LDA), the basic 
difference between Sr 2 VMo06 and Sr 2 VTc06 is in the relative position of the V- and 
Mo(Tc)-states: in Sr 2 VMo0 6 , the V- and Mo-states are separated in energy and form 
two nonoverlapping bands, while in Sr 2 VTc06 these states form one common band 
(Fig. [[]). Both Sr 2 _ :[ ;Y !Z .VMo06 and Sr 2 _ a; Y a ;VTc06 (and other isoelectronic to them 
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Figure 1. Electronic structure near the Fermi level in the local-density approximation 
(total and partial densities of states, and the energy dispersion of the i2g-bands) for 
S^VMoOg (left) and S^VTcOg (right). The position of the Fermi level is shown by 
the dot-dashed line. 

compounds compounds) were intensively discussed in the literature in the context of the 
half-metallic applications. For instance, SrLaVRuOe, SrLaVMoOe, and La 2 VTc06 were 
proposed to be the half-metallic antiferromagnets [HI [161 [E3 HE]. On the other hand, 
the metallic ferrimagnetic (i.e., non spin-compensated) behavior was also suggested for 
SrLaVMoOe [ISl [19]. In general, such conclusions are sensitive to the values of the 
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on-site Coulomb repulsion U, which are used in the calculations. 

The rest of the paper is organized as follows. The method is discussed in Sec. [2j 
Particularly, in Sec. 12.11 we briefly describe construction of the low-energy model for 
the double perovskites and in Sec. 12.21 - the basic idea of the OEP method. Results of 
the model OEP calculations are shown in Sec. EJ Particularly, the obtained electronic 
band structure is discussed in Sec. I3.1[ the screening of on-site interactions due to 
the correlation effects and its consequences on the properties of Sr 2 _ a; Y :r VMo0 6 and 
Sr 2 _ a: Y a ;VTc06 is considered in Sec. 13.21 and the behavior of the total energy and 
stability of the ferrimagnetic states are discussed in Sec. 13.31 Finally, brief summary of 
the work is presented in Sec. HI 

2. Method 

2.1. Construction and parameters of the low-energy model 

The first step of our approach is the construction of the effective low-energy model for 
the t 2 g-states of the double perovskites: 



which is derived in the Wannier basis, by starting from the LDA band structure. The 
Wannier functions are specified by the indices m— xy, yz, or zx, and a stands for the 
spin indices f or \,. The one-electron part t^ m ' of the model Hamiltonian was derived 
by using downfolding method, and the Coulomb (and exchange) interactions V % mm , m „ m ,n 
were obtained by combining the constrained density-functional theory with the random- 
phase approximation (RPA). The method was discussed in the literature, and for details 
the reader is referred to the review article [20]. There are three sets of parameters in 
the model Hamiltonian (pQ): 

(i) the "charge transfer" energies, or the energy splitting between "atomic" ^-levels of 
the V and Mo (Tc) sites, which are described by the site-diagonal part of t^ m . The 
values of these parameters are —391 and 94 meV for the VMo and VTc compounds, 
respectively. Thus, as expected from the LDA band structure (Fig. [1]), for the VMo 
compounds, the atomic V-states are located lower in energy and well separated 
from the Mo-states, while for the VTc compounds, the situation is reversal and the 
splitting is considerably smaller; 

(ii) the transfer integrals, which are described by the off-diagonal elements of t™ m ' with 
respect to the site indices. As expected, the largest parameter is the <i<i7r-transfer 
integral between nearest-neighbor V- and Mo(Tc)-sites, which is 271 (245) meV. 
Another important parameter is the cicia-transfer integral between next nearest 
neighbors of the same atomic type. Its values in the V- and Mo-sublattices of 
the VMo compounds are —72 and —116 meV, respectively. Similar values were 
obtained for the VTc compounds: —73 and —113 meV for the V- and Tc-sublattice, 
respectively; 
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intraatomic Coulomb and exchange interactions in the ^-band, which are screened 
by all other states. At each site of the system, these interactions are fully 
specified by two Kanamori parameters [21]: the intraorbital Coulomb interaction 
U = U mmmm and the exchange interaction J = U mm i m i m for m^m' . In the perfect 
cubic environment, the third Kanamori parameter W = U mmm > m > (the interorbital 
Coulomb interaction) can be related to U and J as W =U — 2 J [20J . Other types 
of matrix elements of ||{7 mm / m >i m "> || are equal to zero. The band-filling dependence 
of U, W, and J , as obtained in the virtual crystal approximation, is explained in 
Fig. [2j Generally, as the number of ^-electrons increases, the effective Coulomb 
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Figure 2. Intraorbital Coulomb interaction U, interorbital Coulomb interaction U', 
and the exchange interaction J for the Sr2_ K Y :E VMo06 (left) and S^-xY^VTcOg 
(right) compounds. 

interactions also increase, reflecting the tendencies of the RPA screening [20]. On 
the other hand, the effective exchange interaction J is less sensitive to the band 
filling. The Coulomb repulsion £Y V is larger for the VMo-compounds, due to 
the specific electronic structure of Sr 2 _a;Y x VMo06 (Fig. [TJ: since the V-band is 
separated and located in the low-energy part of the spectrum, it will be populated 
in the process of the band-filling. It will effectively reduce the number of channel 
available for the RPA screening of W v , for example, due to the transitions from the 
occupied oxygen band to the unoccupied part of the V-band [20]. Thus, U N will be 
less screened in the VMo-systems in comparison with the VTc-ones, where the V- 
and Tc-states form one common band (Fig. [!)). 



2.2. Optimized effective potential method 

In order to solve the model and find the magnetic ground state of the double perovskites, 
we employ the optimized effective potential method. The OEP model works in the spirit 
of the Kohn-Sham density functional theory [9] on the double-perovskite lattice, where 
the kinetic energy 

occ ^ „ 

^Mn = y^y^7; — / dk(c ank \i k \c an]i ) , 

"BZ J 
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and the spin magnetic moment 
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m n 



at the site r are constructed from the one-electron orbitals c CTra k, which are obtained 
from the solution of the auxiliary one-electron problem 

(^k + £>eff) l C cmk) = £ank\c a nk) (2) 

with the potential v° E . In the above notations, is the Fourier image of = \\t mm ||, 
each eigenvector |c CTn k) is specified by its components c^" k in the basis of t 29 -orbitals 
m= xy, yz, or zx of the atomic type r = V or Mo(Tc), n is the band index, and k is 
the position in the first Brillouin zone with the volume fiez- 

The parameters of the potential v° s , which is supposed to be "local" (or diagonal 
with respect to the site indices), are obtained from the minimization of the total energy 

E = E kin + E c + Ex + E corr} (3) 

where the Coulomb {Eq) and exchange {Ex) energies are given by 



and 



err err 

2 / j / j / j mm in in,'" "'mm' "m"m"' 
t era' mm'm"m"' 



Ex 2 ^2 ^2 ^ ; ^mm"'m"m'^mm'^'m"m"'i 

t a mm'm"m"' 



respectively, in terms of the density matrix: 



occ _^ 

*' = J2 7T 



mm' / j f) I crnk crnk' 



jnrfm'i 



The expression for the correlation energy, which is treated in the random phase 
approximation, is given by [221 E3J EM] : 

E COII = J rfqTr jln [l - P(iu, q)tf] [l - UP(iu, q)] + 2P(^, q)tf} , (4) 



where U are the diagonal matrices ||^m'm"m'"ll with respect to the site indices r, 
^ = 1 1 E m 1 m 'm"m"' II i s the polarization in the imaginary frequency: 



pTT > f,'Mflurrr_ / r /k i £<mk ~ £ ™ /k +«J r ™-t 

, "BZ 7 ^ + l^crnk — ^crn'k+qj 



<T n n' 



^ + l^crnk — ^crn'k+qj 



and matrix multiplication in (JJJ) implies the convolution over two orbital indices: 

(UP)mm'm"m"' = zlll' ^mm'U'^ > lVm"m"'- 
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For the cubic systems (and neglecting the orbital polarization effects), the effective 
potential is specified by three parameters: = — A^/2, vj^ = A^/2, v^ ^^ = 
^ - A c M x 0(Tc) /2, and ^ o(Tc) = A ct + A e M x ° (Tc) /2. In the other words, Aj x and A e M x o(Tc) 
is the exchange-correlation splitting between minority- and majority-spin states at the 
sites V and Mo(Tc), respectively, and A ct is the "charge transfer energy" between sites 
V and Mo(Tc). All these parameters are obtained by minimizing the total energy ([3]). 
In practice, first we calculate the total energy on some coarse mesh of these parameters, 
identify the minimum, then continue the calculations around the minimum by using 
finer mesh, and so on until reaching necessary accuracy for the total energy minimum 
(less than 1 meV). 



3. Results and Discussions 



3.1. Electronic structure and arbitrariness of the OEP procedure for half-metallic and 
insulating states 

In the following, we consider two magnetic solutions: ferromagnetic (A^ > and 
^Mo(Tc) ^ ^ ferrimagnetic (A^ > 0, while A^ ^ < 0), and find parameters of 
the effective potential by minimizing the total energy for each case. Total and partial 
densities of states for Sr 2 - x Y 3 ;VMo06, obtained from the solution of the Kohn-Sham 
equations (J2J) with the optimized potential, are shown in Fig. El All FM solutions are 




Figure 3. Total and partial densities of states as obtained in the optimized effective 
potential method for the ferromagnetic (upper panel) and ferrimagnetic (lower panel) 
configurations of Sr 2 _ a; Y :E VMo06. The Fermi level is at zero energy. 
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half-metallic: the Fermi levels crosses the majority-spin band, while the minority-spin 
band is completely empty. In the case of Sr 2 VMo06 and SrYVMoOg, there is a strong 
hybridization between V- and Mo-states, which form the common majority-spin band. 
However, in Y 2 VMo0 6 the majority-spin band is split into the Mo- and V-subbands, 
which are separated by an energy gap. This behavior is related to the large Coulomb 
repulsion LP 7 in Y 2 VMo0 6 (Fig. HJ, which additionally shifts the V-states to the higher- 
energy part of spectrum. Any deviation from the half-metallic behavior and partial 
population of the minority-spin band increases the total energy of Sr 2 _ x Y a .VMo06 and 
thus makes this state unstable. 

The choice of the parameters of the effective potential in the half-metallic regime is 
not unique. For example, any redistribution of states in the unoccupied minority-spin 
band does not affect the ground-state properties. Moreover, any change of A^, A^°, 
and A ct , which does not deform the partially populated majority-spin band (apart from 
the constant shift), does not affect the ground-state properties either. Quantitatively, 
this condition can be formulated as follows: the ground-state properties in the half- 
metallic regime are controlled by the single parameter a = A ct + (A^, — A^°)/2. Any 
change of A^,, A^ , and A ct , which keeps a, does not affect the ground-state properties 
in the half-metallic regime. 

All ferrimagnetic configurations of Sr 2 _ x Y a; VMo06 are metallic. In this case all 
three parameters of the effective potential are uniquely defined. The degree of the spin 
polarization (P) in the ferrimagnetic state is sensitive to the concentration x and the 
definition of P [25] . For example, using the simplest "density of states" definition (see 
Fig. [3]), one obtains P— 0.3, 0.1, and 0.2 for x= 0, 1 and 2, respectively. On the other 
hand, the "Bloch-Bolzmann transport theory" definition yields P— 0.1, 0.1, and 0.5 for 
x= 0, 1 and 2, respectively. The experimental value of P for SrLaVMoOg, measured 
by using the point-contact Andreev reflection technique, is about 0.5 [H]. Nevertheless, 
as we will see below, many of these ferrimagnetic states appear to be unstable and 
calculated spin polarization cannot be directly compared with the experimental data. 

Similar behavior was obtained for Sr 2 VTc06 and SrYVTcOg, where again half- 
metallic FM state competes with the normal metallic ferrimagnetic state (Fig. H]) 

Y 2 VTc06 requires a special attention: it has six electrons in the t 29 -band. The 
system is insulating and all metallic configurations have higher energies. Therefore, 
in the FM state, the parameters A^, A^ c , and A ct can be chosen in an arbitrary 
way. Indeed, any combination of these parameters (provided that the system remains 
insulating) yields the same number of electrons n T = 3, spin magnetic moments 
m T = 3/iB, and the total energy of the system. The correlation energy is exactly equal to 
zero for this particular case. The arbitrariness with the choice of the effective potential 
remains also in the antiferromagnetic state of Y 2 VTcC>6 (Fig. [5]), although the number 
of arbitrary parameters decreases. Indeed, the splitting between V- and Mo-states with 
the same spin affects the ground-state properties. However, the bands with different 
spins can be rigidly shifted relative to each other and, as long as the system remains 
in the insulating state, this shift does not affect the ground-state properties. Thus, 
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Figure 4. Total and partial densities of states as obtained in the optimized effective 
potential method for the ferromagnetic (upper panel) and ferrimagnetic (lower panel) 
configurations of Sr 2 VTc06 and SrYVTcOg. The Fermi level is at zero energy. 
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Figure 5. Total and partial densities of states as obtained in the optimized effective 
potential method for the antiferromagnetic configuration of Y2VTCO6. The Fermi level 
lies in the band gap. 



the properties of the antiferromagnetic Y 2 VTcOg are defined by only two parameters: 
A ct and (A^ — A^), while the third parameter (A^ + A^) can be taken arbitrarily 
Finally, we would like to emphasize that antiferromagnetic Y 2 VTcOe becomes insulating 
due to the large exchange splitting (both at the V- and Tc-sites) obtained in the OEP 
calculations. If the splitting were smaller, the majority-spin V- and Tc-bands would 
eventually merge, yielding the half-metallic antiferromagnetic state, which was obtained 
in the LSDA calculations for the isoelectronic La2VTc06 [IE]. Thus, such half- metallic 
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antiferromagnetic state is probably an artifact of LSDA, and more rigorous treatment 
of the correlation effects will open the band gap in the both spin channels. 

3.2. Screening and the effective Stoner model 

Since the effective potential is local (or diagonal with respect to the site indices) and 
does not affect the orbital degrees of freedom, the OEP scheme has many similarities 
with the Stoner model of magnetism, where the strength of the intraatomic splitting Ag X 
between the minority- and majority-spin states is controlled by the effective parameter 
I T : Ag X = I T m T . Thus, using the intraatomic exchange splitting, obtained from the 
minimization of the total energy, and corresponding values of the spin magnetic moments 
in the ground state (Fig. |6]), one can find the effective Stoner parameter Iqep- All suc h 





Figure 6. Local spin magnetic moments at the V- and Mo(Tc)-sites, as obtained in 
the optimized effective potential method for S^-xY^VMoOg (left) and S^-xY^VTcOg 
(right). 

considerations are valid only for the ferrimagnetic metallic state, where A^ x are uniquely 
defined. The results for Sr 2 _ x Y :r VMo06 are explained in Fig. [7J The effective Stoner 
parameters derived from the OEP approach are of the order of 1.3-1.6 eV at the V-sites 
and 1.0-1.1 eV at the Mo-sites. Similar parameters were obtained for S^-zY^VTcOg: 
7 y EP = 1.27 eV and l£y = 1.06 eV in the case of Sr 2 VTc0 6 , and J V EP = 1.28 eV 
and /qep = 1-18 eV in the case of SrYVTcCV All of them are systematically larger 
that the typical values of I T obtained in the local-spin-density approximation (typically, 
close 1 eV for the 3c?-elements and even smaller for the Ad- and 5<i-elements [26]). 
This example clearly shows the importance of the rigorous treatment of the exchange- 
correlation interactions (beyond LSDA) in the physics of double perovskites. 

Another limiting case is the Hartree-Fock (HF) approximation, which neglects the 
correlation interactions. Then, corresponding Stoner parameter describing the splitting 
in the system of three t2g-orbitals, are given by Jg F = {U T + 2j7" r )/3. They are also 
plotted in Fig. [7J The difference between 7g F and /q EP is the measure of correlation 
interactions and their impact on the magnetic properties of the double perovskites. One 
can see that /q EP is reduced by about 10-20% in comparison with I^ F , which indicates 
the importance of correlation interactions. 

Finally, it is instructive to consider the effect of the static screening alone. The 
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Figure 7. Effective Stoner parameters for S^-xY^VMoOe as derived in the optimized 
effective potential (OEP) method for the ferrimagnetic state, in the Hartree-Fock (HF) 
approximation, and using only static values of Coulomb and exchange interactions in 
the random phase approximation (RPA) . 



corresponding on-site interaction parameters in RPA are given by 

-l A 
U, 

which yields the Stoner parameter /r PA = (W^ PA + 2 l 7 PPA )/3. Such a picture 
can be regarded as the static limit of the GW approximation [271 I2H], where the 
frequency- dependent screened Coulomb interaction £/rpa(^) is replaced by its static 
counterpart Urpa = ^rpa(O) [23] • Typically, it is sufficient to consider only site-diagonal 
contributions to U-rpa, because if U is site-diagonal, the off-diagonal contributions to 
U-rpa are expected to be small [SU] • These results are also shown in Fig. [7J One can 
clearly see that the parameters I PP a are overscreened and strongly underestimated in 
comparison with /qep- This indicates the importance of dynamical screening effects in 
the construction of the optimized effective potential. 

3. 3. Total energy and stability of the ferrimagnetic state 

The behavior of the total energy and of its partial contributions is explained in Fig. El 
For all considered compounds, except S^VMoOe, the ferrimagnetic configuration has 
lower energy. The tendencies towards ferrimagnetism increase with the band-filling. 
The main factor stabilizing the ferrimagnetic state is the gain of the kinetic energy. The 
main factor destabilizing it is the loss of the exchange energy. 

In order to investigate the local stability of the ferrimagnetic phase, we calculate 
the spin stiffness D, which is related to the second derivative of the total energy with 
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Figure 8. Total energy difference between ferrimagnetic and ferromagnetic states and 
partial contributions to this energy difference, calculated for S^-zY^VMoOg (top) 
and Si'2-^Y^VTcOe (bottom). 

respect to the spin-spiral vector q in the ferrimagnetic state: 

E{q) ^ J E(qo) + J D(q-q ) 2 , 

where qo = (vr /a, it /a, it /a) and a is the cubic lattice parameter. More specifically, 
the numerical calculations of E(q), were performed by employing the generalized Bloch 
theorem [3T] and the so-called local-force theorem for the infinitesimal spin rotations. 
The latter states that the total energy change for the small rotations can be replaced by 
the change of the one-electron energies, obtained from the Kohn-Sham equations ([2]) for 
the rigid spin rotations of the effective potential v° s [32]. The theorem is rather general 
and can be proven by considering the rotational invariance of the exchange-correlation 
functions [HJ. Thus, the range of the applicability of this theorem is not limited by the 
local-spin-density approximation, for which it was originally proven, but can be also 
extended to the OEP approach [8], provided that the effective potential ii° s is uniquely 
defined. Particularly, since the choice of v° s is not unique in the case of insulating 
Y 2 VTc06, this strategy cannot be directly applied to this system. 

The behavior of D in the OEP model is explained in Fig. [9j Particularly, D is 
negative in the case of Sr 2 VMo0 6 , Sr 2 VTc0 6 , and SrYVMo0 6 . Thus, the ferrimagnetic 
state in these compounds is unstable with respect to the incommensurate spin-spiral 
state. Since FM state has lower energy in the case of Sr 2 VMo0 6 (Fig. [8]), the spin-spiral 
instability of the ferrimagnetic state may be finally resolved in the favor of the collinear 
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Sr 2 VTcQ 6 SrYVTcQ 6 Y 2 VTc0 6 




Sr 2 VMo0 6 SrYVMo0 6 Y 2 VMo0 6 



Figure 9. Spin stiffness in the ferrimagnetic state calculated using optimized effective 
potentials for S^-xY^VMoOg (filled circles) Sr2_ x Y x VTcOe (empty squares). 



FM arrangement, to make it the true magnetic ground state of S^VMoCV However, 
in Sr 2 VTc06 and SrYVTcOg already ferrimagnetic state has lower energy. Moreover, 
this ferrimagnetic state itself appears to be unstable with respect to the noncollinear 
spin-spiral alignment. Thus, the ground state in this case will be neither ferro- nor 
ferrimagnetic. 



4. Summary and Conclusions 

In summary, using first-principles electronic structure calculations, we have constructed 
the effective low-energy model for the t 2g -states of double perovskites Sr 2 _ a; Y :r VMo06 
and Sr 2 _ x Y :r VTc06. Then, the ground-state properties of the obtained model have 
been investigated on the basis of the optimized effective potential method, by treating 
correlation interactions in the random-phase approximation. The reason for using 
such strategy was twofold. On the one hand, such approach allows one to improve 
the commonly used local-spin-density approximation for the double perovskites and 
incorporate the physics of electron correlations beyond the homogeneous electron gas 
limit. On the other hand, since OEP is formulated as the ground-state method, 
where all the parameters are derived rigorously by minimizing the total energy of 
the system, it gets rid of many ambiguities inherent to the LSDA+£7 method. 
Particularly, although Sr 2 _ x Y a; VMo0 6 , Sr 2 _ a ,Y 2 .VTc0 6 , and other isoelectronic to 
them compounds are frequently proposed to exhibit half-metallic ferromagnetic or 
antiferromagnetic properties, the situation is probably not such optimistic because of 
the electron correlation effects, which affect the intraatomic exchange splitting and 
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makes it substantially larger than in LSDA. Of course, this will revise many predictions 
based on the LSDA. For example, Y 2 VTc0 6 , which was expected to be half-metallic 
antiferromagnet, appears to be insulator. The behavior of other materials is even more 
complicated. For the most of them (except Sr 2 VMoO"6), the half- metallic ferromagnetic 
state is energetically less favorable than the (normal) metallic ferrimagnetic one. 
Nevertheless, this ferrimagnetic state appears to be also unstable with respect to the 
spin-spiral alignment. Thus, the true magnetic ground state of these double perovskites 
will be the spin spiral or some more complicated non-collinear magnetic structure. 

Finally, we have argued that the effective potential (and therefore the one-electron 
band structure) for the half-metallic and insulating ferro- or antiferromagnetic state is 
not uniquely defined. Generally, there is a bunch of parameters, which yield the same 
total energy, occupation numbers and spin magnetic moments at the transition-metal 
sites. 
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